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Abstract 

The quantum electrodynamic theory of the nuclear recoil effect in 
atoms to all orders in aZ is formulated. The nuclear recoil corrections 
for atoms with one and two electrons over closed shells are considered 
in detail. The problem of the composite nuclear structure in the theory 
of the nuclear recoil effect is discussed. 



1 Introduction 



The complete aZ-dependence expressions for the nuclear recoil corrections 
to the energy levels of hydrogenlike atoms were first derived in [1]. These 
expressions consist of three contributions: the Coulomb contribution, the 
one-transverse-photon contribution, and the two-transverse-photon contri- 
bution. For a state a the Coulomb contribution is given by (the relativistic 
units H — c — 1 are used in the paper) 



AE C = AE^+AEW, 
AB« = <«|j£|«>, (1) 
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A ^ 2) = ^/_ oo ^ 5 »(«l[P^c]G(a; + £ a )[p,K]|a), (2) 

where \a) is the unperturbed state of the Dirac electron in the Coulomb field 
of the nucleus, V c = — — is the Coulomb potential of the nucleus, p is the 
momentum operator, S + (u) = -^{oj + ^O)" 1 , G(u>) = {u — H(l — «0)) _1 is 
the relativistic Coulomb Green function, H = ap + (3m + V c . The one- 
transverse-photon contribution is 

AE tr(1) = AESh + AE^ , 

i-(a|(D(0)p + pD(0))|a), (3) 
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r>(uj)G(u + e a ){p,V c })\a), (4) 



where 

D m (u) = -4naZaiD lm (uj) , 

on (I — 1)2,3) are the Dirac matrices, Di m (uu) is the transverse part of the 
photon propagator in the Coulomb gauge. In the coordinate representation 
it is 

D a (»,r) = -±{?^p*M fe + v , Vt («p«H-)-D } . 

47r L r uj z r > 

The two-transverse-photon contribution is 

AE tr{2) = — J_^duj (a\D{u)G{uj + e a )V{u)\a) . (5) 
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An attempt to derive the complete aZ-dependence expressions for the 
nuclear recoil effect was previously undertaken in [2]. Except for the Coulomb 
contribution, the expressions found in [2] are in disagreement with the ones 
given above. A dominant part of this disagreement is caused by technical 
errors made in [2]. If we remove these errors from [2], a discrepancy remains 
in the one-transverse-photon contribution and, in addition, appears in the 
Coulomb contribution. This discrepancy was discussed in detail in [1]. 

Recently, the equations (l)-(5) were rederived in [3,4]. In Ref. [3], it 
was noted that the sum of these expressions can be written in the following 
compact form 

j poo 

AEtot = 2^MJoo dUj{aliP ~ D H) G ( W + ^ p " D H)I°) • ( 6 ) 

The terms AE^ and AE^ can easily be calculated by using the virial 
relations for the Dirac equation [5]. Such a calculation gives [1] 

AE^ s ae® + AE^ = . (7) 

This simple formula contains all the nuclear recoil corrections within the 
(aZ) 4 m 2 /M approximation. The remaining terms (Eqs. (2), (4), and (5)) 
taken to the lowest order in aZ give the Salpeter corrections [6]. Evaluation 
of these terms to all orders in aZ in the range Z = 1 — 100 was done in 
[7]. In particular, it was found in [7] that the complete (in aZ) nuclear 
recoil correction, additional to the Salpeter one, to the Lamb shift (n = 2) 
in hydrogen constitutes -1.32(6) kHz. This value almost coincides with the 
value of the (aZ) 6 m 2 / M correction found in [4,8-10]. 

The complete aZ-dependence expressions for the nuclear recoil correc- 
tions for high Z few-electron atoms were derived in [11]. These formulas 
were used in [7] to calculate the nuclear recoil corrections to all orders in aZ 
for high Z lithiumlike atoms. As it follows from these formulas, within the 
{aZ) A m 2 / M approximation the nuclear recoil corrections can be obtained by 
averaging the operator 

H m = ^E(P* ' P.' - — ("« + ~2 ) ' H ( 8 ) 



s,s' 



with the Dirac wave functions. An independent derivation of this opera- 
tor was done in [12]. The operator (8) was employed in [13] to calculate 
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the (aZ) 4 m 2 /M corrections to the energy levels of two- and three-electron 
multicharged ions. 

In the present paper we generalize the theory of the nuclear recoil effect to 
an arbitrary case of a many-electron atom. In particular, this generalization 
allows one to use as the zero approximation a potential that is different from 
the pure Coulomb field. In addition, it allows one to use the formalism 
in which the closed shells are refered to the vacuum state. In Sec. 2, we 
formulate the basic equations of the method. In Sec. 3 , we apply this 
method to an atom with one electron over closed shells. In Sec. 4, the case 
of an atom with two electrons over closed shells is considered. In Sec. 5, the 
problem of the composite nuclear structure is discussed. 

2 Basic formalism 

Like Refs. [3,4], we will consider the nucleus as a non-relativistic particle. In 
the Schrodinger representation and the Coulomb gauge, the Hamiltonian of 
the whole system is 



where m is the electron mass, M is the nucleus mass, e is the electron charge 
(e < 0), X n is the radius vector of the nucleus, P n = — «Vx„, and fi is the 
magnetic moment of the nucleus. The term —fi ■ H causes the hyperfine 
splitting structure of atomic levels and will not be discussed here. The total 
momentum of the system is given by 



and Pf — J dx(£ t (x) x 7Y(x)) is the electromagnetic field momentum. Since 
the total momentum is an integral of the motion we can restrict our consid- 
eration to the center-of-mass system (P = 0) and, so, can express the nuclear 



H = / rfxV' t (x)[a-(-iV x -eA(x))+/3m]^(x) 




-li ■ H{X n ) , 



(9) 
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momentum in terms of the electron-positron and electromagnetic field mo- 
menta 

P n = -P e - P, = -J dx^ t (x)(-iV x )V(x) - J dx(£(x) x W(x)) . (11) 
Using this equation and the translation invariance we find 
H = J dx^ f (x)[a • (-iV x - eA(x)) + /3m]^(x) 

+ £ / + 5 / *<w + * 2 w> 



i x -y| 

e\e\Z r p e (x) 1 



4tt 7 x 2M 



rfx^ f (x)(-2V x )^(x) 



-y dx(£ t (x) x W(x)) - |e|ZA(0) 



(12) 



Here we have omitted the hyperfine interaction term. The sum of the first 
four terms in (12) is the standard Hamiltonian of the electron-positron field 
interacting with the quantized electromagnetic field and with the classical 
Coulomb field of the nucleus V c = — — (a finite nuclear charge distribu- 
tion can be taken into account by replacing V c with the potential of an 
extended nucleus (see Sec. 5)). The last term in (12) defines the nuclear 
recoil corrections of the first order in m/M. The part of this term containing 
the electromagnetic field momentum ( P/ = /dx(£ t (x) x H(x))) will con- 
tribute only in the first and higher orders in a and, so, will not be discussed 
here. It follows, to the zeroth order in a the nuclear recoil corrections can 
be calculated by adding to the standard Hamiltonian the following term 

Hm = ml ^WHVxMx) / dW f (y)HV y )V(y) 

~ J dx^( x)HVx Mx)A(0) + ^A 2 (0) . (13) 



As is known, for description of an atomic system within QED it is con- 
venient to use the interaction representation in the Furry picture. In such a 
theory the normal ordered form of Hm taken in the interaction representation 
must be added to the interaction Hamiltonian. A simple way to derive formal 
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expressions for the energy level shifts consists in using the technique devel- 
oped in [14,15]. According to this technique, the energy shift AE a = E a —E^ 
of a single level a of an iV-electron atom is given by the formula 

l + &frdEAg a (E) ' ^ 

where the contour T surrounds only the unperturbed level E = E^ (the 
contour is assumed to be traversed counterclockwise), Ag a (E) = g a {E) — 
gf\E),g a {E) is defined by 



/0 

N 



g a (E)5(E-E>) = -—J_Jp\... d p° N dp'l...dp'\ 

x§ (E-p 1 -...-p%)S(E'-p' 1 -...-p'° N ) 

x (u a \G(p'l . . . , p'l-pl . . ,^)T? • • • llK) (15) 

with G is the usual iV-electron Green function 

G(p' 1 ,...,p'° N ;p 1 ,...,p N ) 

= (27c)~ 2N / dx\ . . . dx° N dx\ . . . dx' N 

J — oo 

/ • /0 /0 . . • /0 /O -00 • \ 

x exp (ip ± x ! + ... + ip N x N - ip 1 x 1 - ... - ip N x N ) 
x(0|T^(a;i) . ..^(%gM . . .#n)|0) , (16) 

ip(x) is the electron-positron field operator in the Heisenberg representation, 
ip(x) = ^(x)^f°, u a is the unperturbed atomic wave function, and g^(E) = 
(E — E^) is the zeroth approximation of g a (E). In the first order of the 
perturbation theory we have 

Am Jr 

The Green function G is constructed using the Wick theorem after the transi- 
tion in (16) to the interaction representation. The diagram technique rules for 
G are considered in detail in [15]. Including Hm in the interaction Hamilto- 
nian gives the following additional lines and vertices to the diagram technique 
rules for G. 
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1. Coulomb contribution. 



An additional line (" Coulomb-recoil" line) appears to be 

UJ 



• • 

2-k M -l-oo 



f°° duo 

x y 
This line joins two vertices each of which corresponds to 



where p = — ?V X and k — 1, 2, 3. 
2. One-transverse-photon contribution. 

An additional vertex on an electron line appears to be 




The transverse photon line attached to this vertex (at the point x) is 
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--------- ifZoduDufay). 

x y 

At the point y this line is to be attached to an usual vertex in which 
we have —ie , -foai27r5(uj 1 — uo 2 — u)%) J dy (see [15]), where ot\ (I = 1, 2, 3) 
are the usual Dirac matrices (we note here that in the notations of [15] 
a M = (1, a) and a M = (1, —a) ). 

3. Two-transverse-photon contribution. 

An additional line ("two-transverse-photon-recoil" line) appears to be 

U) 

^^-I-oo duDu(u, x)D lk (u, y) . 

x y 

This line joins usual vertices (see the previous item). 

An important advantage of the approach considered here, in comparison 
with one developed in [1,11], consists in that the present method is suitable 
for arbitrary local potential V(r) (e.g., a local version of the Hartree-Fock 
potential) used as the zero approximation. In addition, the transition to 
the formalism in which the role of the vacuum is played by closed shells 
can simply be realized by changing the sign of iO in the electron propagator 
denominators corresponding to the closed shells. 



3 One electron over closed shells 

Let us consider an atom with one electron over closed shells. In the zero ap- 
proximation the electrons of the atom interact with the potential V(r) which 
can be chosen to include approximately the electron-electron interaction. In 
the formalism with the closed shell states as well as the negative energy states 
refered to the vacuum, the electron propagator is given by 

%x, y) = £ (18) 

where r] n = e n — Ep and Ep is the Fermi energy which is chosen to be higher 
than the one-electron closed shell energies and lower than the energies of 
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the one-electron states over the closed shells. In the simplest case of an 
one-electron atom rj n = e n . 

To find the Coulomb nuclear recoil correction we have to calculate the 
contribution of the diagram shown in Fig. 1. According to the diagrams 
technique rules given in the previous section and [15] we obtain 

(i) _ 1 1 i f°° j (a\Pi\n)(n\pi\a) 



Ag w = -^-i-— / dcoY, ^ ,u,y T7 - (19) 



The formula (17) gives 

M 2rr J-oo ^ oo — e n + ir] n Q 



ae c = J-± r t£ lWM , (20) 



Using the identities 



one can get 



ih = > )+p l (21) 
rho = n5 ^ +v \ < 22) 



' <"IP»-T? E K«IpI»)I 2 - (23) 



2AT 11 ' M , 

£n<£F 

The one-transverse-photon nuclear recoil correction corresponds to the dia- 
grams shown in Fig. 2. A similar calculation gives 

4naZ i f°° ^ ^ [ ( a\p i \n)(n\a k D ik (e a - u)\a) 

lo - e n + ii] n 
| (a|a fc Ajfe(gq - ^)|w)(n|j? < |o) | 



. _, 47raZ i r°° 

J ik (e a -uj)\r 
uj - e n + iry n 



By using the identity 



1 ' 



— e n + «?y n c<j — £ a + iO (a; — e a + iO) (uj — e n + irj n 0) 
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and the equation (21), the expression (24) can easily be transformed to the 
following 

AE tr{1) = AE$ 1} + AE%1 } , 

(2) = W too du)5+{u _ £a) E |Hb^]|n)(n|a fc A fc ( ga -^)|a) 
(a\a k D ik (e a - u)\n){n\\pi,V]\a)\ 



u - e n + ir) n 



(26) 



The two-transverse-photon nuclear recoil correction is defined by the diagram 
shown in Fig. 3. We find 



AE, 



(AiraZ) 2 i 
tr{2) ~ ~1^2^ 



/oo 
-oo „ 



{a\aiDii{e a - Lo)\n)(n\a k Di k {e a - u)\a) 
lo -e n + ir] n 

As it follows from the equations (20), (24), and (27), the sum of all the 
contributions can be written in the following compact form 

1 i f°° 

AE tot = — — / du (a\{pi + 47iaZaiD H (uj)) 
M 271 J-oo 

xG(u + e a )(pi + 4iraZa m D mi (cu))\a) , (28) 

where G(e) = Y,„ — n is the electron Green function. 

In the case of a hydrogenlike atom, the expressions derived here coincide 
with ones given in Sec. 1. 



4 Two electrons over closed shells 

Consider now an atom with two electrons over closed shells (a general case 
of N electrons over closed shells can be considered in the same way). For 
simplicity, we take as the unperturbed wave function the one-determinant 
wave function 

u = ~7n E(- 1 ) P V'Pa(x 1 )Vp 6 (x 2 ) . (29) 
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The nuclear recoil correction is the sum of the one-electron and two-electron 
contributions. Using the diagram technique rules from [15] and the Sec. 
2 and the formula (17) one easily finds that the one-electron contribution 
is equal to the sum of the expressions (28) for a and b states. The two- 
electron contributions correspond to the diagrams shown in Fig. 4-6. The 
two-electron Coulomb contribution is 



/o 

up up 

oo 

1 1 



V p'°-e Pa + iO E-p'°-e Pb + iO 
1 1 

X 

p° - e a + iO E - p° - e b + iO 
x(Pa\p t \a}(Pb\ Pl \b}\. (30) 



Integrating over p°, p'°, and E we get 

AE^ = ^{-l) P {Pa\ Pl \a){Pb\p t \b) . (31) 

A similar calculation of the one-transverse-photon contribution gives 
AE HV = ^E(- 1 ) P [< p a|Pil«><^|a fc D w ( e p 6 - e6 )|6> 

+ {Pa\a k D ki (e Pa - e a )\a)(Pb\ Pi \b)] . (32) 
Finally, for the two-transverse-photon contribution we find 

A< 2 ) } = { ^^U-l) P [(Pa\a k DUep a -e a )\a) 

x{Pb\a m D mi (e Pb - e b )\b) . (33) 

The sum of the two-electron contributions (31)- (33) can be written in the 
following compact form 

AE^f = ^E(- 1 ) P ( Pa l^ + 47r «^A,(^a-e a )|a) 

x (Pb\pi + 4TtaZa m D mi (e Pb - e 6 )|6) . (34) 
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The formulas (31)-(34) coincide with the related expressions found for 
high-Z few-electron atoms in [11] (see also [7]). The difference is only 
the present expressions (31)-(34) are not restricted to the case of the pure 
Coulomb zero approximation. 



5 Composite nuclear structure 

The problem of the composite nuclear structure in the nuclear recoil theory 
was first discussed by Salpeter [6]. In Ref. [6], it was shown that the cal- 
culations based on the assumption that the nucleus is a point Dirac particle 
of electric charge \e\Z and mass M are valid for composite nuclei (indepen- 
dently of the nuclear spin), if the distance between the nuclear levels is large 
compared with the distance between the atomic (electrons + field) levels 
contributing to the nuclear recoil effect. In this section we consider how this 
result can be derived within the approach developed in the present paper. 

Let us assume, for simplicity, that the nucleus is a bound state of a two- 
particle system (e.g., a core with a mass nil and a charge e\ and a valent 
nucleon with a mass m 2 and a charge e 2 ). In this case the sum of the last 
three terms in equation (9) must be replaced by H\ + H 2 + H 3 , where 

Hi = f / rfxp e (x) ( + ), (35) 

47T J \JX — Xi| |X — X 2 |/ 

H 2 = -^( P i-e 1 A(x 1 )) 2 + -^(p 2 -e 2 A(x 2 )) 2 
zmi 2m 2 

+U(x 1 - x 2 ) , (36) 
H 3 = -/!?> ' W(xx) - M 2) • W(x 2 ) • (37) 

Here U (xi — x 2 ) describes the interaction between the nuclear particles (for 
simplicity, we assume that U does not depend of the spins) and /x^ and fi^ 
are the intristic magnetic moments of the nuclear particles. Introducing the 
center-of-nucleus-mass variables 

mixi + m 2 x 2 

An = ; , X„=Xi-X 2 , (38) 

mi + 

we have 

rrii m 2 

Pi = ; Pn + Pn , P2 = ; Pn ~ Pn , (39) 

nil + rri2 nil + m 2 



12 



where P n = -Vx„ and p ra = — V x „ • As in the Sec. 2, we can restrict our 
consideration to the center-of-atom-mass system ( P = P n + P e + P/ = 0). 
So, the total nuclear momentum P n is given by the equation (11). In terms 
of the variables X n and x n the operator H 2 can be represented in the form 

H 2 = H^ + H M) (40) 

where 

H, = ^ + , (41) 

and fj, = (mim 2 ) /(mi + m 2 ). The Hamiltonian if^ describes the intristic 
states of the nucleus. Let us denote the wave function of the nuclear state 
which is under the consideration by a (x n ). The wave function of the whole 
system in the zero approximation is the product of a (x n ) and the atomic 
wave function calculated using the operator (<f) a \Hi\<f) a ) as the interaction 
with the nucleus (we assume here and subsequently that the distance between 
the nuclear energy levels is large compared with the distance between the 
atomic levels). The operator {(j> a \Hi\(f) a ) describes the interaction of electrons 
with the extended nucleus charge. The m/M corrections are calculated by 
perturbation theory. Using the fact that the nuclear size (~ |x n |) is much 
smaller than the atomic size, we expand the vector A in powers of x n . Taking 
into account that a (x„) is of a definite parity we find to the lowest orders 

(4>a\H M \(t>a) 

= jS - ^(Pn • A(X n ) + A(X n ) • P B ) 

I ( e i I e 2 \ A2/V \ 1 ( ei?n2 I e2ini \ 

^\2 mi ~T 2TO2; V n > 2M\ mi ^ m 2 J 

x[p„(x n • V x jA(X n ) + (x n • V x JA(X n )p„]||0 a ) . (42) 

The last term in (42) can be transformed to 

-M^ + ™) J rfx «CK)(x n x p n )0 a (x n )H(X n ) 

= -te<li> + ^^>)W(Xn). ( 43 ) 
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where li and I2 are the orbital moments of the nuclear particles in the center- 
of-nucleus-mass system. Adding this term to the term H 3 gives the total 
operator of the hyperfine interaction — /x • Ti(X. n ), where /x = /x^ + /xj, 1 ) + 
nf^ + /x( 2 ) is the total magnetic moment of the nucleus and /x^ = ^-lj. Due 
to the operator H M contains the term 

^ = -(— -— )pn-A(X n ), (44) 

there is a contribution of the order m/M from the second order of the per- 
turbation theory. For a state a of the whole system we have 

AE , = t gi e 2 x2 ( a |p ra ■ A(X ra )|w)(w|p w ■ A(X w )|a) 

Assuming that the energy difference between the nuclear state which is under 
the consideration and the other nuclear states contributing to the sum in the 
equation (45) is large compared with the corresponding energy differences 
between the atomic (electrons + field) states which give a dominant contri- 
bution to AE' a , we replace E a — E n in Eq. (45) with e a — e n , where e a and e n 
are the nuclear energies. Using the identity p n = i^H^ x n ] we find 

v mi m 2 ' 2 n ^a ea ~ €n 



x[(0 a |[iJ M ,X„]|0„)(0 n |p n |0 a ) + (0 a |Pn|0n)(0n|[^,X n ]|0 a )] 
= (— --) 2 ($a|A 2 (X n )|$ a )^V 

v mi m 2 y 2 



X[(0 a |x n |0 n )(0 n |p n |0 a ) - (4>a\Pn\4>n)(4>n\Xn\4>a)} 



$ a |A 2 (X n )|$ a )-/x(0 a |[x„,p n ]|(/. a ) 



v mi m 2 y 2 

= 4(— )V(^|A 2 (X„)|$ a ), (46) 
2 v mi m 2 7 

where |0 a ) is the nuclear wave function and |<& ) is the atomic (electrons + 
field) wave function. Combaining this term with the equation (42) we find 
that the nuclear recoil correction of the first order in m/M is defined by the 
operator 

Hm = ||-^^(Pn-A(X n )+A(X n )-P n ) 
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(e 1 + e 2 ) 2 2 
+ 2M A (Xn) 

1 (P n -|e|ZA(X n )) 2 . (47) 



2M 



So, within the approximations made in the derivation of this result, the effects 
of the composite nuclear structure can be neglected. 
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Figure 1: Coulomb nuclear recoil diagram. 



o 

a b 

Figure 2: One-transverse-photon nuclear recoil diagrams. 
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Figure 3: Two-transverse-photon nuclear recoil diagram. 
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Figure 4: Two-electron Coulomb nuclear recoil diagram. 



19 



Figure 5: Two-electron one-transverse-photon nuclear recoil diagrams. 



Figure 6: Two-electron two-transverse-photon nuclear recoil diagram. 
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